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Abstract. We study the entanglement between the 2D vibrational motion and 
two ground state hyperfine levels of a trapped ion, Under particular conditions this 
entanglement depends on the parity of the total initial vibrational quanta. We study 
the robustness of this quantum coherence effect with respect to the presence of non- 
dissipative sources of decoherence, and of an imperfect initial state preparation. 
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1. Introduction 



Over the last few years we have witnessed a very rapid development of cooling and 
trapping techniques both for neutral atoms and for ions |jl|, |[ ||, |, ||, 0. These 
progresses have made it possible to perform sophisticated experiments wherein several 
examples of couplings between a few bosonic and fermionic dynamical variables have 
been realised p|, ^, |l^, [Tl|, |12|. If, in fact, an ion confined in a miniaturised Paul trap is 
exposed to suitably configured laser beams, the 3D harmonic motion of the ionic center 
of mass gets entangled with the internal degrees of freedom. Some peculiar aspects of the 
vibronic response stemming from such a dissipation-free situation have been successfully 
exploited for realizing experimentally Fock states, coherent states, squeeezed states and 
Schrodinger cat-like states 10, 11], presenting theoretical schemes for engineering 
several nonclassical state |1^, |14|, |15|, ^ [1^, |19|, [2^, 0, |23[ , implementing quantum 
logic gates |2^, realising tomographic reconstructions of the density matrix of the 



29|, and, more in general, discovering quantum effects characterising 



system 

the ionic quantized oscillatory motion j^O], |3^, |3^. Many hamiltonian models 

have been reported so far in the literature to describe physical properties of trapped 
ions. Some models explore physical situations wherein only the ionic motion along a 
specified direction of the trap is effectively influenced by the presence of the laser beams. 
When the induced vibronic coupling is instead extended over two (three) independent 
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directions of the trap, then one refers to the correspondent physical scenario as to a 2D 
(3D) trapped ion. 

In this paper we study the motion of an ion isotropically confined in the radial 
plane of a Paul microtrap when it is irradiated by a properly chosen configuration of 
external laser beams. We show that there exist experimentally interesting conditions 
under which the dynamics of this system may be exactly treated [^. Exploiting this 
fact we find that the center of mass vibrational motion and the electronic degrees of 
freedom of the ion becomes entangled, in a way which is very sensitive to the parity 
of the initial number of vibrational quanta. Then we study the robustness of this 
quantum effect with respect to various experimental imperfections. In fact, even if 
dissipation of the vibrational motion is completely negligible, non-dissipative sources 
of decoherence associated with fluctuations of external parameters may be important. 
We study these effects, and those of laser intensity fluctuations in particular, using the 



model-independent formalism of Ref. ^j. We find that, using standard values for the 
strength of the laser intensity fluctuations, the parity-dependent entanglement effect is 
not washed out by non-dissipative decoherence. The effect of an imperfect initial state 
preparation of the vibrational degree of freedom is also studied. 



2. Parity dependent entanglement effect 
2.1. Derivation of the effective Hamiltonian model 

Consider a two-level ion of mass M confined in a bidimensional isotropic harmonic 
potential characterised by the trap frequency v. Let's denote by Kuoq the energy 
separation between the ionic excited state |-|-) and its ground state |— ) and assume 
ujq 3> v. Indicate by a (a^) and h {f)^) the annihilation (creation) operators of vibrational 
quanta relative to the ionic center of mass oscillatory motion along the x and y axes of 
the radial plane of the trap respectively. 

We assume that the ion is driven by two 7r-out of phase laser beams, applied 
along the two orthogonal directions x and y with an angle of 7r/4 relative to the x 
and y axis respectively. The two lasers have equal frequency ujl = ujq — 2z/, intensity 
Eq and wavevector modulus k^. After performing the dipole and the rotating wave 
approximations, the hamiltonian model of the system can be cast in the form 

H = Tiu{d)d + h^h) + + [dt^~^ (x, y, t)a_ + h.c] (1) 

where d is the appropriate dipole matrix element, assumed real for simplicity, and 
e^~\x,y,t) = £'Q[e*(^^*~^^^) _ g«('^Lt-fcLj/)j represents the negative frequency contribution 
to the resultant electric field acting upon the atom. The two-level atomic degrees of 
freedom are described in terms of the inversion operator = |+)(+| — |~)( — I and of 
the two transition operators a-t = |±)(=f|. In a reference frame rotating at the laser 
frequency ujL,the hamiltonian H assumes the form = Hq + Hi with 

Ho = hu{d^d + b^b) + —a, (2) 
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Hi = /ifi{[e-'''(^+^') - e-^^^^+^'^l^T^ + h.c.} (3) 

where 



5 = u;,-u;, = 2u, fi=^, V = kL]J^ (4) 

are the atom-laser detuning, the Rabi frequency and the so-called Lamb-Dicke parameter 
respectively. We have indicated with A {A'' ) and B {B'^ ) the annihilation (creation) 
operators of vibrational quanta along the directions x and y respectively, defined by 

i = -L(a + i) B = ^^{h-a) (5) 
In the interaction picture, the interaction Hamiltonian assumes the form 

- f; {-ly+^^^^^B^ (fit)'e-*('=-i-2)^_ + h.c] (6) 

where the resonance condition 5 = 2v has been used. Under the assumption v ^ Vt, 
we may perform the rotating wave approximation neglecting in Eq. (P) all the terms 
oscillating at multiples of the harmonic trap frequency v. Finally, using Eqs. (|^) it 
is not difficult to show that, in the Lamb-Dicke limit r/ <^ 1, the effective interaction 
Hamiltonian in the interaction picture can be written as [|l^, ^ 



Hint = % (ab (T+ + h.c}j (7) 
with g = Vtrfe'''^^/'^. 

2.2. Experimental set up 

For realizing the Hamiltonian model we make use of a three level electronic system 
as shown in Fig. 1. Such a system is similar to the one currently used in the experiments 
performed at NIST [10|. The weak transition |— ) ^ |-|-) is simultaneously driven by 



two Raman laser configurations with different detunings (Ai and A2) from the electronic 
level |3) and with frequency differences ujl = i^Li — = = — ^^l^ = "^o ~ 2i/. 
For sufficiently large detunings, such that the intermediate state |3) can be adiabatically 
eliminated, only the electronic levels |-|-) and |— ) are involved in the dynamics. In 
this case the two-photon Raman transitions between the states of interest are formally 
equivalent to narrow single photon transitions [1^]. The Raman laser fields are chosen 



to have wave vector differences k l = k — k and k = k — k pointing in 
the X and y directions respectively. Under these conditions, in the resolved sideband 
regime and in the Lamb-Dicke limit, the dynamics of the system is described by the 
Hamiltonian (0) provided that we make the following substitution g = i7/j?7^e~'^«''^ with 
Qr = VLl^VLl^/Z^i = Ql'^Ql'J A2 and r^/j = ^Jhkl/2Mv = ^Jhkl,/2Miy. The dependence 
of the Lamb-Dicke parameter on the modulus of the wave vector difference allows to 
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Figure 1. Three- level electronic system of the trapped ion and scheme of Raman 
couplings 



vary r^/j by changing the Raman beams configurations. In this way it is possible to 
obtain very small values of the Lamb-Dicke parameter as required in our scheme. 

For the scope of the paper it is important to underline that measurements of the 
electronic states are experimentally realizable with very high efficiency by means of 
the so called quantum jumps tecnhique ||, |10|, [ill. The method essentially consists 
in coupling the ground electronic level |— ) to a third auxiliary level \r) by means of 
another laser beam. The transition |— ) ^ |r) is chosen to be a dipole allowed one. 
In this conditions, the presence of fiuorescence detects the ion in the electronic ground 
level and its absence in the excited level. With this method state detection efficiency 
approaches 100%, even if a modest number of scattered photons are detected. 



2.3. Dynamics 

Let's denote with \na, rib) = \na)\nb) the product of the vibrational Fock states along 
the X and y directions and suppose that the initial state of the ion has the form 



|vl/(0)) = |r = l,jo = f)|-) = ^E, ^ 

k=0 \ 



\N~k,k)\-) 
^J2P,\N-k,k)\-) (8) 



N 



k=0 



The vibrational state |r = 1, jo = y) belongs to the class of the so called SU(2) coherent 
states defined as 



1 y^2i 



/ r'\2j-k,k) 



(9) 
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where r G C, 2j G N. The states \N — k,k) appearing in Eq. (||) are eigenstates of the 
operator (ji'^d + b'^bj all pertaining to the eigenvalue N = 2jo representing the initial 
total number of vibrational quanta. In particular, as pointed out by Gou and Knight 
the SU(2) coherent state corresponding to the specific value t = 1, |r = l,jo) of 
a bidimensionally confined ion, coincides with the vibrational Fock state with = 2jo 
quanta along the direction x forming an angle 7r/4 relative to the x axis. Since number 
states for the ID motion of an ion in a rf-trap have been already experimentally realized 
(see |ll|), it is evident that the initial state of Eq. (||) could be easily prepared 
experiment ally. 

Then, if at t = the laser fields realizing the coupling between the 2D vibrational 
motion and the internal two-level system described by the Hamiltonian model of Eq. (J^) 
is turned on, at any subsequent time instant t, the state of the system can be written 

as 



with 



and 



\^it)) = \v-m-) + \Mm+) (10) 

N 

= Y.Pk cos(/fet)|iV - k, k) (11) 

N-1 

= -^T.Pk Mfkt)\N -k-l,k-l) (12) 

k=l 

where 

fk = 2g^{N-k)k (13) 

are the Rabi frequencies. Eq. (p]0| ) shows that, starting from the factorized initial 
state l^'(O)) of Eq. (^, the Hamiltonian model (|^) generates entanglement between the 
external and internal degrees of freedom of the trapped ion, giving rise to interesting 
dynamical consequences. In order to appreciate the meaning of this assertion, we focus 
our attention on the time evolution of the vibrational entropy 

S^it) = -Tr[p,{t) In p,(t)] , (14) 

Py being the reduced density operator describing the 2D vibrational motion of the ion. 
A straightforward calculation gives 

S,{t) = - In [c(t)^W [1 - c(t)]^-^Wl (15) 



where 



N 1 



c{t) = J2\Pkfcos\fkt) 2 



N 



Y,\Pk\cos{2fkt) 

k=0 



(16) 



fc=0 

Exploiting an analytical method based on the analysis of the Rabi frequencies, it has 
been recently demonstrated in Ref. that, starting from a total vibrational excitation 



number A^, there exist an A^-dependent time instant at which the internal and external 
degrees of freedom of the trapped ion are disentangled (c(t) = 1 or c(t) = 0) or maximally 
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Figure 2. Time evolution of the probability to find the ion in the internal ground 
state P-igt) for iV = 9 (a) and 10 (b) respectively. 

entangled (c(t) = 1/2). More in detail, it has been analytically proved that if ^ 1 
is odd, at the time instant to = one has c(to) = 1 (c(to) = 0) if (A^ — l)/2 is 

even (odd). This implies that Sy{io) reaches its absolute minimum, indicating that the 
internal and external degrees of freedom manifest a marked tendency to disentangle 
each other. On the contrary, if ^ 1 is even, at the time instant ig = one has 
c{te) = 1/2 and therefore the vibrational entropy {tg) reaches its maximum value. 
This means that the vibrational and electronic degrees of freedom are, at this time 
instant, maximally entangled. Fig. 2 displays the time evolution of the probability 
P-{t) = c{t) of finding the ion in its ground state (which can be measured with the 
technique described at the end of subsection 2.2) according to Eq. ([16|); Fig. 3 instead 
shows the vibrational entropy ^^(t), as given by Eq. (p^. In both figures, the cases of 
an initial total vibrational number = 9 (a) and A^ = 10 (b) have been considered. 

These figures illustrate in particular the existence of a A^-dependent time instant 
at which the system under scrutiny exhibits different quantum behaviours dependent 
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(a) 




Figure 3. Time evolution of the vibrational entropy Sv{gt) ior N = 9 (a) and = 10 
(b) respectively. 

on the parity of N. The interesting physical aspect is given by the pecuhar nonclassical 
sensitivity to the granularity of the initial total number of vibrational quanta N. The 
physical origin of this intrinsically quantum behaviour stems from the specific two- 
boson coupling mechanism envisaged in this paper. In the rest of the paper we shall 
analyse the robustness of this quantum behaviour with respect to various experimental 
imperfections as the fluctuations of the intensity of the coupling lasers and a non-ideal 
initial state preparation 

3. Influence of non-dissipative decoherence sources 

In the previous section we have assumed perfect unitary evolution and initial state 
preparation for the system. This assumption is partially justified by the fact that a very 
good control of the quantum dynamics has been now achieved for trapped ions and also 
because they are well isolated from their environment. In fact, many experiments have 
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shown that dissipative effects on the vibrational motion can be considered completely 
negligible for long intervals of time. Moreover, spontaneous emission is also practically 
irrelevant because the two internal states are two hyperfine levels of the ground state. 

However, decoherence effects have been nonetheless observed in the motion of a laser 
driven trapped ion In fact, even when the entanglement with the environment 

is negligible, fluctuations of some classical parameter of the system may cause non- 
dissipative, phase-destroying effects. In the trapped ion case, phase decoherence is 
mainly caused by the fluctuations of the Rabi frequency, which are induced by the 
laser intensity fluctuations, making the coupling constant g fluctuate. A quantitative 
explanation of the decoherence observed in the Rabi oscillations of Ref. P], in terms of 



a fluctuating laser pulse area, has been recently provided by Bonifacio et al. in Ref. |37 
using a model-independent formalism able to describe non-dissipative decoherence 
phenomena due to the fluctuations of classical parameters or internal variables of a 
system. 

In the model-independent approach of Ref. the dynamical quantities in the 
presence of non-dissipative decoherence can be obtained simply through an average 
over an appropriate probability distribution of the fluctuating parameter. In the 
trapped ion case, the random parameter is the positive dimensionless random variable 
^(^) — lod^diO^ which is proportional to the laser pulse area. In fact, the quantity 
g{^) is a stochastic coupling constant, whose fluctuations have to be traced back to the 
laser intensity fluctuations. The time evolution of a generic dynamical quantity 0{t) 
becomes therefore the averaged quantity 

0{t) = / dAP{t,A)0{A) . (17) 
Jo 

The probability distribution P{t, A) is the Gamma distribution function 

e-A/ar (A/^r) , . 



Pit, A) 



gr T{t/T) 



obtained in Ref. |33] by imposing the semigroup condition for the averaged density 
matrix of the system. The physical meaning of the parameters g and r in Eq. (^) 
can be easily understood by considering the mean and the variance of the probability 
distribution (Hi 



{A) = gt (19) 
a^A) = {A') - (A)' = gHr, (20) 

implying that g has now to be meant as a mean coupling constant, and that r quantifies 
the strength of A fluctuations. In Ref. [0 the decay of the Rabi oscillations observed 
in Ref. is well fitted assuming r ~ 1.5 ■ 10^*^ sec (see also |0). In this case, the 
relevant experimental timescale t is much larger than r and in this limit the Gamma 
distribution function can be well approximated by the following Gaussian distribution 

P^t,A) ^ —l—expL^-^--^] . (21) 
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Figure 4. The probability difference AP^ defined in the text, at the time instant 
t, as a function of the laser pulse area fluctuation strength parameter r (expressed in 
seconds), in the case of a total initial vibrational quanta = 9 and for g — 10^ Hz. 



The experimentally observed quantity is the probability of finding the ion in its ground 
electronic state -P-(t) = c(t). In the presence of the fluctuating coupling constant g{^) 
caused by the intensity fluctuations of the coupling lasers, this probability has therefore 
to be averaged over the Gaussian distribution (^). Using Eq. ([TBI) we flnd 



where the quantum coherent oscillations at frequencies fk are now exponentially 
decaying because of the non-dissipative decoherence. 

To quantify the robustness of the parity-dependent entanglement with respect to 
the various experimental imperfections, we have considered the quantity AP^ , deflned 
as the difference between the values assumed by P-{gt) in correspondence to (odd) 
and N+1 total initial vibrational excitations, at the "intermediate" time i = With 
no imperfections and iV > 1, one has i ^ ie ^ to, and P^='"^'^{gi) ~ 1, P^+^igt) ~ 1/2, 
so that AP^ ~ 1/2 whatever is. We emphasize that such a result highlights the 
fact that the probability P-(t) of flnding the ion in its electronic ground level at t = t 
depends on being even or odd. As long as the effect of phase decoherence increases, 
the difference between the odd A^ case and the even A^ case tends to zero, because in 
both cases the system state becomes the same statistical mixture and AP^ tends to 
zero. 

In Fig. 4 we show the dependence of the probability difference AP^ as a function 
of the laser pulse area fluctuation strength r, in the case when A^ = 9 and g = lO^Hz. 
One can see the expected decay to zero of AP^ for increasing r (i.e., increasing laser 
intensity fluctuations); in particular one has a well visible transition from a quantum 
behaviour at small r to a completely decoherent behaviour at larger r, taking place at 
a threshold value Tf/j ~ 3 • 10~^ sec. Since current experimental situations correspond to 
T < 10^* sec, this means that the parity-dependent entanglement effect described in the 




(22) 
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preceding section can be experimentally detected even with the unavoidable presence of 
non-dissipative decoherence sources, as for example laser intensity fluctuations. 



4. Imprecision in the initial state preparation 

Another important condition for the observability of the parity effect reported in this 
paper is the possibility of preparing the initial vibrational Fock state |A^) along the 
direction X' at tt/4 radians with respect to X. It is however not difficult to guess that 
at the end of any scheme aimed at preparing a vibrational ionic Fock state containing 
exactly vibrational quanta, the state effectively reached is a vibrational density matrix 
Pv{0) containing also diagonal terms diffeent from | A^) (A^| in the total vibrational energy 
eigenbasis. Of course, such a mixture of initial odd and even total vibrational quanta 
counters the possibility of observing the parity dependent entanglement. However, 
we shall show that the degree of preparation efficiency required in order to keep the 
parity effect in the dynamics of our system appears realistically in the grasp of the 
experimentalists. In other words we shall prove that the occurrence of parity effects, 
although attenuated, turns out to be compatible with the initial state preparation fidelity 
presently achievable for ionic vibrational states. 

We assume that at t = the initial vibrational state can be described as follows: 

oo 

p.(0) = E Are-(™-^)'/2^V,0)(m,0|, (23) 

m=0 

J\f being the normalization constant. Equation (^) is a gaussian weighted mixture of 
pure Fock states |m, 0)(m, 0| along the direction X', centered on a prefixed number A^ 
of vibrational quanta and controlled by a width parameter A. It is immediate to see 
that when A ^ Pv{0) tends to |A^, 0) (A^, 0|. The form given to p^(0) suggests the 
introduction of the parameter 

^ = 1 _ ^ = 1 _ e-V2A^ (24) 
Pn 

describing the efficiency of initial state preparation, p^+i and appearing in Eq. (|^) 
are the probabilities of having A^ + 1 or A^ quanta in the initial vibrational Fock state. 
With the help of equations (^) and (|23|) , it is not difficult to verify that, in the 
presence of both laser intensity fluctuations and imperfect initial state preparation, 
the probability to find the ion in the internal ground state |— ), P_{gt), assumes the 
form 

1 



P^gt) 



1 + E Xe-^rn-N)y2A- ^ ^-89^r(m-k)kt cos(4y'(m - k)kgt) 
m=0 A;=0 



(25) 



As we have done in the preceding section, we quantify the robustness of the parity- 
dependent effect in terms of the probability difference AP^ defined above. In Fig. 5 
AP^ is plotted as a function of the preparation efficiency rj for three different values of 
the fluctuation strength parameter, r = 10~^ sec (solid line), r = 10~^ sec (dotted line) 
and T = 10~^ sec (dashed line), again in the case N = 9, g = 10^ Hz. As expected, one 
has a steady increase of AP^ for increasing preparation efficiency. In particular, if we 
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Figure 5. The probability difference AP^ at the time instant i as a function of the 
efficiency of the initial state preparation ij (see Eq. (|2^)) for r = 10~^ sec (solid line), 
r = 10~^ sec (dotted line) and r = 10"'' sec (dashed line). The other parameters are 
as in Fig. 3. 

assume the realistic parameters rj = 0.9 and r = 10^^ sec, one can see that the parity 
dependent entanglement effect is still visible, although attenuated by almost 40%. 

5. Conclusions 

In this paper we have analyzed how some experimental imprecisions influence a 
nonclassical effect we have brought to light and discussed in a previous paper P^ , 
namely the parity dependent entanglement effect. According to this effect by changing 
only one vibrational quantum in an initial condition characterized by a total number of 
oscillatory quanta ^ 1, the state of the system, at approximately the same instant 
of time te ^ to, drastically changes. Indeed at this time instant the vibrational and 
electronic degrees of freedom of the ion are maximally entangled (i.e. maximum reduced 
entropy) for even N and disentangled for odd A^. 

The parity dependent quantum effect has a nonclassical origin since it stems from 
the granularity of the oscillatory energy of the center of mass of the ion. For this reason 
we believe it's worth being experimentally verified. We have thus studied the robustness 
of the effect with respect to various experimental imperfection. Since single trapped ions 
are essentially dissipationless systems, we have considered the effect of non-dissipative 
decoherence sources as the intensity fluctuations of the coupling lasers, and also the 
effect of an imperfect initial state preparation. The main result of the paper is that the 
parity-dependent entanglement is quite robust against these experimental imperfections 
and that it can be experimentally seen using presently available technology. 
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